arXiv:1508.06964v2 [hep-ph] 22 Dec 2015 


Electron in three-dimensional momentum space 


Alessandro Bacchetta, 1, 2 ’[2] Luca Mantovani, 1, 2, [j] and Barbara Pasquini 1, 2 >|f| 

1 Dipartimento di Fisica, Universita degli Studi di Pavia, 1-27100 Pavia, Italy 
2 Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 
(Dated: December 23, 2015) 

We study the electron as a system composed of an electron and a photon, using lowest-order per¬ 
turbation theory. We derive the leading-twist transverse-momentum-dependent distribution func¬ 
tions for both the electron and photon in the dressed electron, thereby offering a three-dimensional 
description of the dressed electron in momentum space. To obtain the distribution functions, we 
apply both the formalism of the light-front wave function overlap representation and the diagram¬ 
matic approach. We perform the calculations both in light-cone gauge and Feynman gauge, and we 
present a detailed discussion of the role of the Wilson lines to obtain gauge-independent results. We 
provide numerical results and plots for many of the computed distributions. 

PACS numbers: 14.60.Cd, 12.20.-m, 13.60.Hb 


I. INTRODUCTION 

The electron is an elementary field of quantum electrodynamics (QED) and, as such, it is a structureless object. 
When we study the electron “structure”, what we are in fact doing is probing the quantum fluctuations of the theory. 
Due to the Heisenberg uncertainty principle, the electron can fluctuate into a virtual electron-photon pair, carrying 
the same quantum number as the electron, i.e. e —» ey —> e. The virtual photons can themselves break up into pairs 
of virtual electrons and positrons and, as a result, an isolated electron is in reality surrounded by a composite virtual 
cloud of photons, electrons, and positrons in all possible combinations, consistently with the quantum number of the 
electron. Therefore, the “bare” electron effectively becomes a “dressed” electron. Photons, electrons, and positrons 
can be interpreted in this case as partons contained in the original electron, in analogy to the partonic structure of 
hadrons. 

If one of the virtual states of the electron cloud interacts with a probe, the parton content of the electron is re¬ 
solved and the electron reveals its structure. Such behavior is an important aspect of quantum field theories, and 
similar phenomena arise in many different situations. For example, in QCD, the quark splitting to quarks and gluons 
drives the scaling violations in hadronic structure functions pj. Analogously, the QED radiative corrections to the 
electron can be recast in the formalism of structure functions (see, e.g., Refs. mm)- Other possible descriptions of 
the parton structure of the electron can be obtained in analogy to the QCD formalism of parton correlation functions, 
which can be parametrized in terms of different types of distribution functions (see, e.g., Refs. HE]) known as gen¬ 
eralized transverse-momentum-dependent parton distributions (GTMDs), generalized parton distributions (GPDs), 
transverse-momentum-dependent parton distributions (TMDs) and collinear parton distribution functions (PDFs). 

These distribution functions allow us to view the dressed electron from a new perspective. For instance, it is possible 
to introduce the notion of the shape of the electron 130 , to analyze how the spin of the electron is split into the 
spins and orbital angular momenta (OAM) of the partons in the electron’s cloud fOHTTI . and to compute all possible 
correlations between spin, momentum and position of the partons. Moreover, the study of QED parton distribution 
functions can help to shed light on several formal aspects of their QCD analogues. 

In this work, we will focus on the study of TMDs, which describe the distribution of QED partons in the dressed elec¬ 
tron as a function of their longitudinal and transverse momentum, and therefore give access to the three-dimensional 
picture of the dressed electron in momentum space. At leading twist there are eight TMDs, for both the active electron 
and photon, characterizing the strength of different spin-spin and spin-orbit correlations of the partons and electron 
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target. Three of them survive when integrated over the transverse momentum, giving rise to the familiar collinear 
PDFs. Formally, the transverse-momentum dependence in the distribution functions arises from the definition in 
terms of bilocal operator matrix elements off the light cone. To comply with gauge invariance, one has to introduce 
nontrivial gauge links, which compensate for the nonlocality of the fields p2HT6] . The Wilson line accounts for the 
(infinitely many, in principle) photons that can be absorbed or emitted by the initial or final state, and it is the source 
of interesting phenomena like single-spin asymmetries mm- 

In a diagrammatic approach, the gauge-link contribution can be studied via a perturbative expansion in the coupling 
constant. Thanks to its Abelian nature, QED is simpler to deal with compared to QCD, and can therefore serve as a 
useful testing ground for full perturbative calculations. In this article, we will work in the leading-order approximation, 
corresponding to the exchange of a single gauge boson between the active parton and the spectators. In particular, we 
will use QED Feynman rules for the eikonal propagator, and we will explicitly show the equivalence of the resulting 
TMDs in different gauges, namely Feynman and light-cone gauge. 

Besides the diagrammatic approach, a useful framework to study parton distributions is the representation in terms 
of overlap of light-front wave functions (LFWFs) [19]. The LFWFs of the electron corresponding to the minimal Fock 
state component of an electron-photon pair are computable explicitly from time-ordered perturbation theory in the 
light-cone gauge. Several applications using the electron LFWFs have been discussed, focusing in particular on the 
study of the distributions in impact-parameter space as obtained from the GPDs [31 113 IMIS]- Here we will extend 
the study with LFWFs to the TMDs, of both the electron and photon partons. 

The diagrammatic approach and the LFWF overlap representation are used as alternative, but equivalent, descrip¬ 
tions to emphasize different aspects of the physical content of the TMDs. The diagrammatic approach allows us to 
unravel more explicitly the role of gauge invariance in the definition of the TMDs, while the representation in terms of 
LFWFs, being eigenstates of the parton light-front helicity and OAM in the light-cone gauge, exposes more explicitly 
the information about the spin-spin and spin-orbit correlations encoded in the TMDs. 

The outline of this work is as follows: In Sec. [il] we set up the general framework for the TMDs of the electron, 
introducing the definition of the correlator function and discussing the role of the Wilson line in QED. In Sec. |III| we 
derive the leading-twist TMDs using the diagrammatic approach in the Feynman gauge, at leading order in a. In 
Sec. |IV| we sketch the derivation of the LFWFs for the electron-photon Fock state component and discuss the LFWF 
overlap representation of the TMDs, showing the equivalence of the results within the diagrammatic approach. In 
Sec. [V] we study in some detail the role of the transverse components of the Wilson line in the light-cone gauge, using 
different prescriptions for the regularization of the light-cone singularity from the photon propagator. In Sec. VI we 


turn our attention to the photon contribution to the dressed-electron TMDs. In Sec. Em we present the results for 
the electron and photon TMDs, discussing the role of the spin and parton OAM in shaping the transverse-momentum 
dependence of the different TMDs. Finally, we conclude with a section summarizing our findings. 


II. GENERAL FRAMEWORK 


Let us consider a dressed electron of mass m, spin S and four-momentum P, composed of a bare electron of 
momentum k and a photon of momentum P — k. If we write a generic vector in light-cone coordinates as 
a M = (a + , a - , aj_), we set the momentum of the inner electron as k = {xP + , fc - , kj_), with the longitudinal momentum 
fraction x defined as x : = k + /P + . We begin with the definitions related to the distribution functions of the internal 
electron; those for the photon will be discussed in Sec. [VT] In analogy to the QCD case, if we probe the internal 
structure of the electron via the scattering of a virtual photon (see Fig. [TJ , the most general expression for the 
correlator function between initial and final electron states is given by [U [2LI1T281 


: = / (§i elk< 5 l^(°M(0l p, s ) ■ (i) 

In Eq. (JT]) ip{x) is the electron field, while the indices i and j refer to Dirac space and will be omitted from now on. 
If we integrate 0 over k , we obtain the transverse-momentum-dependent correlator, which enters the definition of 
the TMDs: 
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FIG. 1: Handbag diagram for the dressed electron in QED at order a (left), separated in its upper and lower part (right), the latter being 
the diagrammatic representation of the electron-electron correlator <f>. 


If we further integrate ([ 2 ]) over kj_ , we reduce ourselves to a correlator function that depends only on the longitudinal 
momentum fraction x, giving access to the collinear PDFs: 

$(x-,P,S) = J dk-d 2 k ± $(k;P,S) = j ^e ifc -«<P,^(0)V>(£)I^S> 

From a diagrammatic point of view, the correlator corresponds to the lower part of the handbag diagram in Fig. [T] 
In order for the bilocal products of field operators in the correlator functions to be gauge invariant, we have to 
insert a Wilson line between the two electron fields [HUS 12312S1 ED]. Indeed, the QED Lagrangian 1 

~2qed = i>(i$ - m)ip - ^F‘ JV F IIV + e^ 7 ^ipA^, (4) 

where the photon field tensor is — d v A lx {^), is invariant under the following local gauge transfor¬ 

mation of the fields: 

^( 0 +^( 0 . ( 5 ) 

However, this transformation does not leave the correlator $ invariant, since 

^(0)V>(£) -t Vi(0)V>(0 e - ie [ a(0) - a( « )] . (6) 



The standard way to ensure gauge invariance of the correlator $ is to insert in our definition a gauge-link operator 
Lf(o,£) connecting the fields evaluated at different points. Let us define, in general ITB) . 


%l, 6) :=ex P 



(7) 


where the integral is meant to run along any path connecting to £ 2 - Acting with the gauge transformation ([5]), the 
Wilson line will turn into 
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which is exactly the transformation that we need to render $ gauge invariant. 

Since different processes require different paths along which the Wilson line must be evaluated [HI ;18] (29l |3Q], 
we are now going to introduce the specific path that enters the definition of the transverse-momentum-dependent 
correlator |14l . Consider a generic vector v^ = (v + ,v~,v±) in light-cone coordinates; for our purposes we can keep 
the plus component i> + fixed at 0, since the fields in ([ 2 ]) are taken along the light cone at £ + = 0. We move in 


1 Here, and in the following, we take e > 0, meaning that the electron charge is equal to — e. 
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the v -v_l plane and define Wilson lines along the v and v_l directions as follows: in order to connect the points 
(0, a~, cj_) and (0, b~, c_l) along the minus direction, we introduce the longitudinal gauge link: 


U v (a ,b ;c±): = Vexp 


-ie / dr/ A + ( 0,7? ,cj_) , 

J a~ 


(9) 


while in order to connect the points (0,c , aj_) and (0, c , bj_), along the transverse direction, we use the transverse 
gauge link: 


W ± (a_ L ,b J _;c ) 


V exp 


r b_L 

ie / dr)± ■ Aj_(0, c~, r) ± ) 
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( 10 ) 


In order to link the points 0 = (0 , 0j_) and £ = (£ , £j_) appearing in our correlator, we follow a path along the v 
direction all the way to infinity, then get to £_l along the transverse direction and finally go back to , as shown in 
Fig. m To this purpose, we define 

o,<x>)'- = U v (0~, oo _ ; 0j_)^/ J_ (0_ L , oo_l; oo _ ) , (11) 

W(oo,0 ■■ = U ± {oo±,€±;oo~)U v (oo~,£“;£_l) . ( 12 ) 


We therefore define the gauge-invariant transverse-momentum-dependent correlator as: 
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It is clear that if we choose the light-cone gauge A + = 0, the longitudinal gauge link U v reduces to the identity; in this 
gauge, only the transverse gauge link survives, and this is crucial to maintain the gauge independence under residual 
gauge transformations, which can be fixed by imposing appropriate boundary conditions 1 1 31 . In Feynman gauge, 
instead, the longitudinal gauge link is nontrivial; however, the transverse gauge link does not give any contribution, 
because the boundary condition for the transverse component of the photon field is Aj_(oo) = 0. In the following, we 
will discuss the differences between the two choices of the gauge in more detail. 

The insertion of a gauge link brings about some changes to the correlator from a diagrammatic point of view mm- 
Wilson lines physically correspond to final- and initial-state interactions, i.e. they account for the (infinitely many, in 
principle) photons that can be absorbed or emitted by the incoming or outgoing particles. Some examples of this kind 
of interaction are shown in Fig. [3] where we gather all the a-order contributions to the correlator (diagrams (a) to 
(f)), along with a case of a higher-order contribution (diagram (g)). Diagrams (a) and (b) are actually of zeroth order 
in the gauge link, as they refer, respectively, to the handbag diagram we already discussed and to the correction due 
to the self-energy of the bare electron; diagrams (c) and (d) contain one gauge photon and they are hence described 
using the first-order approximation of the Wilson line in the coupling constant e, while (e) and (f) are of second order 
in the expansion of the gauge link. We also notice that, in diagram (d), the photon loop is responsible for the virtual 
vertex correction, which contributes only at the end point x = 1, kj^ = 0j_; this is true also for the self-energy in 
diagram (b). Diagrams (e) and (f) are instead proportional to n 2 _ [with nt = (0,1, OjJ] and can be neglected, unless 
the direction of the gauge link is not taken along the light cone. 

We remark here that a proper definition of TMDs must deal also with the end points and the occurrence of infrared 
and rapidity divergences, which can be regularized, e.g., by taking n_ to be off the light cone. These are critical 
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FIG. 2: Wilson lines connecting two points at 0 = (0 . Oj_) and £ = (£ ,£_l). 
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FIG. 3: Contributions to the correlator at different orders in a and in the gauge link. Notice that also the Hermitean conjugates of all 
diagrams should be taken into account. 


issues in proving factorization theorems and introducing well-defined TMDs, both for QCD and for QED. Careful 
investigations of these issues can be found in several papers and books (see, e.g., [iI5 | [5TH57] ). 

The correlator d>(fc;P, S) in 0 can be written as a combination of elements of a basis in Dirac space T> : = 
{l, 7 ^, 7 5 , 7 5 7 Let us consider from now on the case of the transverse-momentum-dependent correlator 
$(x, kj_; P, S) defined in ([2]) [2511551155] . We can choose a reference frame where the incoming dressed electron has 
no transverse momentum, i.e. P = (P + ,P~, 0_l), with P _ = m 2 /2P + for the on-shell condition. Moreover, since it 
must be P • S = 0, S' 2 = —1, we can set S = (S Z P + /m, — S z m/2P + 1 Sj_) and therefore replace the dependence on the 
covariant four-vector S with the dependence on the three-vector S : = (Sj_, S z ), with S 2 = 1. 

Introducing also the light-like vector n+ = (1,0, Ox) and the antisymmetric tensor, 


= -e 2 , 1 = 1 


= e 2 ± = 0 , 


(14) 


one finds, in the leading-twist approximation, 
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where i,j = 1,2 and o([ = (0,0, aj_). In the above equation we introduced the leading-twist TMDs of an electron 
inside the dressed electron, /f, fif, gf L , gf T , h\ Tl h>i£, hi£, hi e , that depend on (a:, kj_). As far as the notation 
is concerned, the letters /, g and h refer to unpolarized, longitudinally polarized and transversely polarized internal 
electrons; L (T) refers to the spin of the dressed electron being along the longitudinal (transverse) direction; the “_L” 
symbols signal an explicit dependence on transverse momenta with an uncontracted index. Out of the eight TMDs 
in Eq. (15 1, the so-called Boer-Mulders function hi e m and the Sivers function m are T-odd, i.e. they change 
sign under naive time reversal, which is defined as usual time reversal but without interchange of initial and final 
states, while the other six are T-even. 

Let us take the trace of the correlator multiplied by an appropriate Dirac operator and define 


$[ r l(x,k ± ;S):=iTr[<I>(x,k ± ;S)r] , 

with T € T>. It turns out that, at the leading twist in 1 /Q, these Dirac space projections are given by 
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FIG. 4: Diagrams contributing to the TMDs in Feynman gauge, for x yf 1 and k x Ox- 


where 


h\ : = h e 1T 


2 m 2 


nr 


( 20 ) 


The different TMDs in Eqs. (|T7|)-(|T9j) can be isolated by taking proper spin configurations of the dressed electron, 
while the Dirac structure T in the field correlator projects on different polarization states of the internal bare electron. 


III. CALCULATION OF TMDS IN FEYNMAN GAUGE 


In this section, we derive explicit analytic expressions for the TMDs. We stress that, from now on, we will exclude 
the end point x = 1 , kx = Ox from our analysis, since it would require a proper regularization procedure which is 
beyond the scope of the present work. This implies that we will neglect in our analysis all self-energy and virtual 
vertex corrections discussed in section |H| (diagrams (b) and (d) in Fig. [3]). 

An important statement can be made straightforwardly about the T-odd TMDs and h^ e . Had we tried to 
derive them starting from the correlator without the gauge link, i.e. the one in Eq. (|2|, we would have found that 
they are vanishing. This, eventually, can be explained by the fact that only diagrams containing a loop at one side 
of the cut (such as diagrams (d) to (g) in Fig. [3]) can potentially give a non zero contributions to T-odd TMDs M- 
The only possibility to have such a loop is either to consider diagrams of higher perturbative order, or consider the 
end point x = 1, kx = Ox- Therefore, already at this stage we are allowed to conclude that at order a and for x 1, 
kx 7 ^ Ox, we have: 

hi e (x,k i) = 0, (21) 

frr( x > kx) = 0 ■ (22) 

As we mentioned, from a diagrammatic point of view, the correlator ([l]) corresponds to the diagram shown in the 
right-hand side of Fig. |T] it is possible, then, to evaluate it with usual Feynman rules and take its trace with a proper 
matrix F, in order to recover the analytic expressions of the TMDs. Introducing for greater clarity the Dirac indices, 
the Feynman rules give for diagram (a) in Fig. [ 4 ] 

\M\ 2 = ((P - k) 2 ) 5 {{k + q) 2 - m 2 ) , (23) 

where we factor out the transition matrix corresponding to the lower part of the diagram, 


:=e 2 d^(P-fc)«s,fc(P) (7, 


k 2 — m 2 — it 


kj 


+ m 


—7m) u s ,i{P) , 

ie J il 


(24) 


from the one corresponding to the upper part, where the interaction between the parton electron and the virtual 
photon takes place, 


M v p : = e 2 (/l(q)(^ + ^ +■ 

We can henceforth write for the correlator, setting the Dirac indices according to (|T]) 

mS) ij = j£y i 6(( p -k) 2 )M'° w . (25) 
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FIG. 5: Handbag diagram with a photon loop separated into its upper and lower part, with the latter corresponding to the correlator. 
The two factorized subdiagrams are meant to be “attached” by joining the crosses with the crosses and the circles with circles. 


The sum over photon polarization vectors d^ v takes different forms depending on the gauge. We choose the Feynman 
gauge first, and hence use 

(PK(P) =-<T- (26) 

A 


As an example, we would like to show the result for the TMD fi(x, kj_), which is the distribution for unpolarized 
electrons in an unpolarized dressed electron. From Eq. since it is fyf = 0, we see that /f can be recovered by 
taking the trace of <f> 7 + (with a factor 1/2) and including also an average over the spin state of the target electron. 
We rewrite also: 


6{(P~k) 2 ) = 


2(P — k)+ 


6[k~ -P~ + 


ki 


2 (P - k)+ 


where, from the on-shell condition, P = to 2 / (2P + ) (working in a frame of reference where Pi = 0j_)- Hence, it is 
easy to find that the contribution of diagram (a) of Fig. [4] to the TMD /f is 


/i°(*,kl) = X! / dk ~ 


8(2t r) 3 (P- fc)+ 


5 [ k — P + 


ki 


(-e 2 ) 


2(P — k)+ J (k 2 — m 2 ) 2 


(2t r) 3 


k^_ + m 2 (l — 4x + x 2 ) 


(fi 2 (x, k 


a). 


[u s { P)7m(^ + TO ) 7 + (# + m^usCP)] 


(27) 


where we drop the imaginary part of the electron propagator since the internal electron cannot be on-shell, and we 
define the function ip{x,\s? L ) as 


<d*,ki) 


e \J x(l — x) 

[(1 — x) 2 m 2 + ki] 


(28) 


This is not actually the final result, since we have not yet considered the longitudinal gauge link, which is not 
vanishing in the Feynman gauge and brings about two more diagrams to be evaluated, namely diagram (b) and (c) 
of Fig. [4j Applying the Feynman rules to diagram (b), one finds 2 


\M\ 2 = -e A S (l 2 ) 




k 2 — m 2 — ie 




! T jf + / + m 
(k + q + l) 2 — m 2 + ie 


/ A (qbs(P) . (29) 


In order to separate the term corresponding to the lower part, this time we need to apply the so-called eikonal 
approximation BUGS], which consists of considering as relevant only the ” component of the momentum of the 
fermion after the interaction with the virtual photon; moreover, it is assumed that the emission or the absorption 
of a “soft” photon, with small momentum l, does not alter the fermion momentum. Under this approximation the 
propagator of the struck fermion becomes 


$ T (J[ T / T tyi 
( k + q + l) 2 — m 2 + ie 


21 + 


( 30 ) 


2 For convenience, we set the momentum of the dressed electron as P = k + l. 
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where we use 


[k -f- q -(- Z) (A) -f- q Z) 4- , (A) T q + Z) A)j_j T (^_Li T Z_l i > (A; + (7 H- Z) m , (k -\- q) Z , 

along with the on-shell conditions (fc + g) 2 = m 2 , Z 2 = 0 . In this way we are able to factor out a lower and an upper 
part in the amplitude, since 


{ft + f + m)7 P 7 + ~ (A; + q) 7+7^7 +S p = 2{k + q) 7+ ~ 2{ft + + m) 


(31) 


where in the first and in the last step we applied again the eikonal approximation. The modified propagator described 
in Eq. (30) is shown in the right-hand side of Fig. |s} the eikonal approximation is taken into account by using the 
modified Feynman rules for the eikonal propagator and vertex described in Fig. [6] 


Z 

1 

1+ + it’ 



= —ien 




FIG. 6: Feynman rules for the eikonal propagator (left) and vertex (right), when appearing on the left-hand side of cut diagrams. Notice 
that the momentum in the propagator is equal to the momentum of the photon attached to vertex. 


The resulting transition matrix for the lower part of the diagram in Fig. [5] is now 

ft + m 


M low = - 

t3 1+ + it 


-g^n-v-UsA p ) 7i 


k 2 — m 2 — it 


«s,i(P) ■ 


(32) 


kj 


Let us focus our attention on the denominators first. As we already noticed, the imaginary part of the electron 
propagator can be dropped, and we are left with 


1 
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l + + it k 2 — m 2 l + (k 2 — m 2 ) 


— in- 




(33) 


where the term with S ( l + ) will be canceled by the contribution of diagram (c) of Fig. [4J since the imaginary unit 
from the photon vertex changes sign, being on the right-hand side of the cut. Therefore, we are allowed to drop the 
it also in the photon propagator. If we derive the correlator corresponding to the amplitude (32), as we did in (25) 

the con 

fi( x i kj_) = 


starting from (24), we obtain that the contribution of diagram (b) to the TMD /f is 

k]_ + m 2 (l — x) 2 " 1 


(2tt) 3 


( I -*) 2 


tp (x,ki) . 


(34) 


A completely equivalent result is found by computing the contribution f-f to /f from diagram (c); henceforth we 
have fi = fi, and we can finally write as the final result 


fi(x, kf) = A“(x, ki) + fi(x, ki) + /f(x, ki) = 


(2tt) 3 


(1 + x 2 )k5_ + to 2 (1 — x) 4 
x(l — x) 2 


if 2 {x, k^_) . 


(35) 


The same procedure can be applied to all the TMDs; as a result, we find that the leading-twist TMDs, at (x, kj_) A 
(1,0), are 
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hiT(x,k ± ) = 0 , 

ht(x,kl) = 


ki 


(27 r ) 3 (1 — x ) 2 
h~i e {x, ki) = fir(x,k\_) = 0 . 


<p 2 (x,k 2 ± ) = 


2 e 2 


xk 2 ± 


(2 tt) 3 (1 — x) [kj_ ± m 2 ( 1 — a;) 2 ] - * 


(40) 

(41) 

(42) 


Note that the results for the T-even TMDs in Eqs. (36)-(41) correspond to the leading-order predictions in the 
quark-target model of Ref. [32] , up to the QCD color factor. The Abelian nature of QED shows up in the results for 
the T-odd TMDs, which are vanishing at variance with the quark-target model. 


IV. CALCULATION OF TMDS IN LIGHT-CONE GAUGE 


A. Overlap representation 


We now switch to the light-cone gauge, aiming to recover the same results that we found in the Feynman gauge. We 
will also apply a different approach, adopting the formalism of light-cone quantization and thus obtaining an overlap 
representation of TMDs in terms of LFWFs mmm- 

We already remarked that the longitudinal gauge link is simply an identity in the light-cone gauge; we also neglect, 
at least at this level, the contribution coming from the transverse gauge link, that we will study separately later. 
Therefore, the expression we should start with is ([2]). Taking the trace of a generic Dirac operator T is equivalent to 
rewriting [25] 


$I r l(z,k ± ;S) 




£+=0 


(43) 


We expand the eigenstate |P, S) of the electron in the light-cone Fock space H3] truncating the expansion to the 
two-particle Fock state of a photon and an electron. Since we are going to use light-cone quantization, we need 
to decompose the dressed electron states in terms of light-front lrelicity states |P, A). The relation between the 
eigenstates of canonical spin S and the eigenstates of light-front helicity A can be found, for instance, in Ref. 23- In 
this representation, we then have 


|P, A) = Z\e) + (ey) , 


(44) 


where A = ± stands for the values of light-front helicity ±1/2. In the above equation Z is a renormalization constant 
that we may simply take to be 1 for our purposes (it would become relevant only when considering the end point 
x = 1, kx = 0_l). The one-particle state is simply: 

|e) = 6l(P)|0) , (45) 

with b]y (P) the creation operator of an electron with light-front helicity A and momentum P = (P + ,Pj_). 

In order to find the exact expression of the multiparticle Fock states, one should solve the eigenvalue equation for the 
full QED Hamiltonian. However, by considering the interaction potential small with respect to the free Hamiltonian, 
one is allowed to apply time-ordered perturbation theory and, for the two-particle state at first order in the coupling 
constant e, one finds [i¥ 


|ey) 


dxd 2 k± 

2(2tt) 3 ^/x(1 — x) 


^A.A^ukOley; p e ,A; p 7 ,A 7 ) . 

a,a 7 


(46) 


Here A = ±, A 7 = ±1 stand for the internal electron light-front helicity ±1/2 and photon polarization, respectively, 
and p e , p 7 refer to the momenta of the partons 3 : 

p e = (zP + ,k_L) , p 7 = ((1 - P)P + ,-kj_) . (47) 


3 One should pay attention not confusing the \e r y) state, which describes the dressed electron as a 2-body system in Eq. (j44j, with the 
state of two free quanta |e 7 ; p e , A; p 7 , A 7 ) . 
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The wave functions 'I ; A A (x,kj_) appearing in (46) are given by 


^VwObk-O = <^(a;,ki) u x (p e )^£l m (p 7 )«a(P) 


(48) 

with the function y?(x, kj_) defined in Eq. ( [28| . In the light-cone gauge A + = 0, they are the LFWFs for the |ey) 
component of the dressed electron. For a dressed electron with light-front helicity A = +, one has explicitly 4 jjSj 


'F+ + i(x,k L ) = \ - 


2 ( k X iky) 


x (1 — x) 


<p( x i k 


±) > 


_r(a:,kj.) = ~y/2x 


(h x + iky') 


</?(>, k l) , 


(1-x) 

'Ft !+ i(x,k J _) = ^m(l - x) ip{x,k 2 ± ) 


*± _i(*,k x ) = 0 . 


Assuming instead A = —, one has 






kj_) = 0 , 


i(aJ, kj_) = \ -m( 1 — x) (p(x, kj_) 




,+i 


(x, kj_) = V2x^-f- - l ^y- <p(x, k5_) , 

(1 -x) 


k ±) = -\l~ (l _ %) 


2 (fox H - 'iky) 




Notice that the LFWFs in Eqs. (49)-d56k are related by 


*A,A>>kx) = (-l)CA-^)/ 2 +^4/ : ^: A (x,kx). 


(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 


If we use Eq. (|44|) in the correlator (|43|), and take the proper projections for each F matrix, we can obtain the 

rhich reads 

*A',V ( x , (p e )ru A (p e )4'A,A 7 (x, k -0 • 


LFWF overlap representation of the correlator, which reads 

1 


$[ r l(x,k ± ;A) = 


4(2t t) 3 xP+ 


(58) 




An alternative, yet equivalent, way to derive the TMDs explicitly is to represent the correlator in the basis where 
we consider the light-front helicity of both the dressed and the internal electron, and we treat them symmetrically. 
We can therefore define the light-front helicity amplitudes [45l : 


1 


$A',A';A,A(z, k _L) := (P, A' \b{, (x, k ± )5 A (x, kj_) \P, A) 


(59) 


with 


N : = [2(2tt) 3 ] 2 x<5 (3 ) (0) . 


Inserting in (|59|) the Fock-state expansion p4|), we obtain 

1 

m 

A 


$A',A';A,A(a;, k x) = 3^3 ( a > k -L)^A,A 7 (x, k -L) ■ 


( 60 ) 


4 There is a sign difference with respect to Ref. [9], due to a different choice of the sign of the charge e. 
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The light-front helicity amplitudes are parametrized by the following combinations of TMDs 


$ = 


/ 

Ufi+9 e 1L ) 



+ Pit) 

K \ 


t( ih ± e + h i£) 

|(/f ~ 9\l) 

k t 7 | p 

2m7 

t^(*/ir — Pit) 


- Pit) 

k R U±e 

2m 2 ,1 1T 

s(/i e - 

Pit) 

-t^i e + hi£) 

V 

K 

+ 9it) 


- Kl) 

\ift + g e 1L ), ) 


(61) 


where k^L : = k x ± ik y . The row entries are (A', A') = (++), (H—),(—b),(-), while the column entries are 

(A, A) = (++), (-1—), (—h), (—)• 

Using Eq. ( [58] ) and the projections or equivalently Eq. ( [60| along with ( [61] ) , one obtains the following 

results for the LFWF overlap representation of the contribution of the internal electron to the T-even TMDs: 

fl(x,k 2 ± ) = I^A^^ik-L)! 2 = [l^+,+i(a: ) k_ L )| 2 + |^+_ 1 (a;,k_ L )| 2 + |^+ +1 (a;,k_ L )| 2 ] , (62) 


PilOu^) = 

gl T (x,k 2 ± ) = 

hi£(x,k 2 j_) = 


1 


2(2tt) 3 

m 


4(27r) 3 k2_ ^ 


^lT^.kj.) = 


2(27r) 3 k2_ 

m 2 

2(27r) 3 k-y 


[|d't, + i(a;,k i )| 2 + |^_!(a;,kj .)| 2 - |^± i+ 1 (a:,kj_)| 2 ] , 

J2 [ k R * I,a 7 (*> kO'F- (x, kj_) + k L ^ (*, k ± )$+ ^ (x, k ± ) 


[ kR (*’ k -0*+.X, (*> k ±) + k L 'F+^ (x, kj_)tf+ A-r (x, kj.) 


E[ 


h e i(x,k 2 ± ) = 


2(2t r) 3 


*i,+i 


k R ^-,A T (*> kj_)^ +iA ( x , k J_) + k L tf +iA (x, kx)^_ )A (x, kj_) 


{x, k_ L )\[ , _^ +1 (a;, kj_) + ^(z, k ± )T_ _ 1 (x, kj.) 


(63) 

(64) 

(65) 

( 66 ) 

(67) 


The analytic results found by inserting the explicit expressions (49)-(56) of the LFWFs indeed coincide with the 


results (|36|)-(|42j) obtained in Feynman gauge and using Feynman diagrams. 

The LFWFs are eigenstates of the total spin of the partons in the z direction, J z = A/2 + A 7 , and of the total 
OAM L z = A/2 — J z , as required by angular momentum conservation. For the two-body LFWFs of the electron, one 
can have only L z = 0, and L z = ±1. They are commonly labeled as S and P waves, respectively, although this is an 
abuse of language as partial waves should refer to L and not L z . The LFWF overlap representation in Eqs. (62)-(67) 
allows one to disclose the different contributions to the TMDs from the spin and OAM configurations of the partons. 
In particular, /-f, pf L and (if are all diagonal in the OAM, while the remaining TMDs contain different interference 
terms between S and P waves, and therefore involve a transfer of OAM between the initial and final electron states. 
The interplay between the different partial waves in the TMDs will be discussed in more detail in Sec. m 


B. Feynman diagram approach 


It is interesting to follow the diagrammatic approach in light-cone gauge also, in order to explicitly discuss the 
gauge invariance of TMDs. In this case, since the longitudinal gauge link is trivial and we are neglecting for the 
moment the transverse gauge link, we only need to evaluate diagram (a) of Fig. [4] the amplitude is still given by (231, 
but this time we take the following sum over the photon polarization vectors 5 : 


d^( P ) ■■= (pk( P ) = -<r + — (pV/ + P v) 


( 68 ) 


If we factorize the amplitude separating the upper and lower parts, we actually get three contributions for the lower 
part of the diagram: 

Ad'ow = A + B + C , 


5 Notice that the third term of the polarization sum, proportional to p 2 , is omitted since in our case the photon is on shell. 
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with 


A : = -e us (P)j, 


B : = 


^(Pb -1 


ft + m 


(. P-k)+ 

f + m 


ft + TO 


-jK)«s(P) , 


-7 + «s(P) • 


(69) 

(70) 

(71) 


It is now easy to recognize that the term A is equivalent to the contribution of Eq. (24) evaluated in the Feynman 
gauge. Let us focus on B ; using the Dirac equation {f — m)its(P) = 0, we can rewrite 

,2 + m 2 ~ ^) U s(P) = ,2 + m 2 ( m - ft) U s(.P) = -«s(P) ■ 

Hence we find that term B , which comes from the extra term in the polarization sum in light-cone gauge, actually 


coincides with the contribution of Eq. (32), treated with the eikonal approximation. The same applies also to term 


C 1 which in turn corresponds to diagram (c) of Fig. [4] in Feynman gauge after eikonal approximation. 

We can now also easily check that the term of the correlator that is related to diagram (b) and (c) of Fig. [4] are 
actually vanishing in the light-cone gauge. If we still use the eikonal approximation, but substitute the sum (|68[) 


instead of (261 for the polarization vectors in (32), the latter becomes: 


M low = -- 


-Mfi(P)7^ 


(. P-k)+ + ie 
which contains the Dirac structure, 

^7 + 


-r + 


i 


(p-k)+ 


((P - kyn v _ + (P - k) v nt) 


+ m 


k 2 — to 2 — ie 


n_ M u s (P) , (72) 


1 


(P - k)+ 


7 /ib) (P — k) v n- v + 


(P-k)+ 


— ft)n 2 _ = 0 . 


Within the diagrammatic approach in the light-cone gauge we recover the same results for the TMDs we obtained in 
the Feynman gauge derivation of Sec. |III| 


V. TRANSVERSE GAUGE LINK 


In the previous sections we performed the derivation of the TMDs for the electron first in Feynman gauge including 
the contribution from the longitudinal gauge link, then in light-cone gauge but neglecting the Wilson line; as a result, 
we found perfect agreement between the two gauges, as it should be, since the TMDs are gauge-invariant objects. 
However, one may wonder if there is a contribution of the transverse gauge link. We actually neglected a correction 
to our results that is illustrative to analyze, although it will show up only at the end point x = 1. 

As it was noticed by Ji et a 1. USES], the transverse gauge link does have an important role in taking into account the 
final-state interactions when working in light-cone gauge. The point is that the four-momentum of the intermediate 
photon must be integrated over, according to cut-diagram rules. This means that, when we use the sum over the 


polarization vectors (68), we actually should regularize the singularity due to the p + factor in the denominator. We can 


in fact do this in different ways, because we have the freedom of choosing a certain prescription for the regularization 
of the denominator; the following choices are possible: 


Retarded : 


Advanced : 


1 


p + + ie 




n+ — 


Principal value : —— 

p+ 


1 


p + + ie p 


+ - 


ie 


(73) 

(74) 

(75) 


It is obvious that our final result must be independent of the choice of a specific prescription, and it must also be 
the same we found in Feynman gauge, where the problem of regularizing the denominators does not arise and the 
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transverse gauge link does not contribute. Nonetheless, if we come back to our evaluation of the TMD /f (x, kj_) from 
Feynman diagrams, we see that it is not the case if we make the replacement (73)-(75) in (68) and then put the latter 
in the amplitude (24). 

Let us first take the retarded prescription as an example: term A in (69) would be unchanged, but the denominator 
of terms B and C in (70) and would be modified according to (|73|), giving 


B = 


(P — k)+ + ie 
and similarly for C. We rewrite, as usual: 

1 


Ms(P)7 _l 


Z + m 


Z + to 


k 2 — to 2 — ie k 2 — m 2 + ie 


’-Z)ua(P) , 


(P - k)+ + ie 


= P.V. 


(<pA yf) - - Q+) ■ 


(76) 


(77) 


This leads us to the final result for the TMD (including all three terms): 


/i ret (*,ki) = 


1 


k 2 t 


(1 + x 2 ) 


(27t) 3 [ x(l — x ) 2 

„ 2 i 


= fi(x,k ± ) - i- 


2 {l-xY 

TO - 

X 

<5(1 — x) . 


ip 2 (x,ki) - 




———|-to 2 (1 — x) ip 2 (x, kj) 5(1 — x) 

1 — X 


Wki 

Similarly, we obtain for the advanced prescription: 


f? v (x,ki) = m x ,ki) + i- 


i 


<5(1 — x) . 


Wk 2 

In the principal-value prescription, instead, the denominators of B and C remain unchanged because: 


(78) 


(79) 


1 1 

-;-1-r-■ 

p + + ie p + — ie 


P.V. 


1 


— i7T(5(p + ) + P.V. 


1 


+ iTrS(p + ) 


1 


and hence the final result coincides with the one found in Feynman gauge: 

fr(x,k 2 ± ) = mx,ki). 


(80) 


We expect all results to be the same when adding also the transverse gauge-link contribution. As we will show, the 
extra terms proportional to <5(1 — x) in the retarded and advanced prescriptions are compensated by the contribution 
of the transverse gauge link, while the latter is vanishing in the principal-value prescription. We remark that this is not 
what one finds usually in the case of the nucleon (see e.g. [12]), since the retarded prescription is the one corresponding 
to having (oo) = 0 and hence the transverse gauge link should vanish in that prescription. However, we stress 
that in the present discussion we have not taken into account all the corrections to the distributions at the end point 
x = 1, and this might explain the difference with the existing calculations in QCD. 

We now turn to the evaluation of the transverse gauge link in the LFWF overlap representation. We start again 
with the integrated version of the correlator, with the states in the light-front helicity basis: 


®(x, kj_;A) : = 


ik< 


(27 


e lk< (P,A^(0)W ( 0 iOO )W (oOi £)^(0|P,A) 


«+=o 


(81) 


We express the Wilson line in the light-cone gauge, using (10) approximated according to 

W( 0 ,oo)^(oo,£) - 1 + ie [ dr)± ■ Aj_(0, r/±, oo) = t + ie f dtA±(0,t^xp 
J Oi Jo 


3 ) -£l 


(82) 


We neglect the identity in (82) and we focus on the gauge-link terms. If we switch again to light-cone quantization 
and replace the dressed electron state \P, A) according to (44), in principle we will get four different combinations 
of states to consider. However, the term containing the gauge link between two jey) states is vanishing, because the 
photon coming from the final- (or initial-) state interaction cannot couple to the spectator photon. The only relevant 
contributions will then come from the terms where the matrix element is between the states \e) and ley); in this 
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situation, the gauge photon which interacts with the bare particle in the initial (or final) state coincides with the 
photon in the dressed electron. Therefore, we are left with the following two terms, contributing to the correlator: 


I:= jJ ^(2^)^ W(°) eft Aj_(0,Cl,oo) • £_l ) V’(O)le'y) 

d(~d 2 $, ± ik< 


II :=ie 


(2tt) 3 


e M e 7hK0) / dt A_l(0,^_l,oo) ■ W’(O)le) 


«+=o 


«+=o 


(83) 


(84) 


Let us now restrict to the case of the TMD /®; we will begin with the evaluation of I. After some manipulations and 
resolving the integral over t , it can be split into two terms, I = I a + 7b, with 


Ia: = ~ 


e f d£ d 2 £± 


2(27r) 3 J 2(27r) 3 ^x'(l — x') 


A _gjfe-Sg-iPe-Jg-ip+oo 


x H ~ fP '!r ^ »A(P)7 + ^A(p e )4'A ;A7 (a: / ,k'x) 

P7j_ * S_L 


a,a,a 7 

e f dt~d 2 €_L 
+B : = x 


£+=o 


(85) 


efa/^k' 


- 1 g*fe-?g-*Pe'Cg-jp+OO e i$j_-p 7 


2(2tt) 3 J 2(2 7 r) 3 v V(l - x') 
x J2 ±X ^ P7 l ^ u A (P)7 + UA(Pe)^iA 7 (x',k' ± ) 

A,A,A 7 P7J.A-L 


£+=o 


where p e = (x'P + , k' ± ) and p 7 = ((1 — x')P + , —k' ± ). The sum appearing in both I a and Ib is 

T: = u A (P)7 + u A (p e )u A (p e )j M u A (P)£^ (p 7 )e±^(p 7 )- = -4P + p e± ■ £ ± (l - 2^") , 

A,A,A, ' P 1 ' 


( 86 ) 


(87) 


where we used the explicit expression of the LFWFs in terms of the spinors and polarization vectors. Equation (85) 
consequently becomes 


I A = -2 eP+ 


d£, d 2 £ l f dx’d 2 k'x 

2(2tt)3 J 2(2tt) 3 v V( 1 - x ') 


0 ik-t e -i Pe -Z _ e -ioop+ +2x 'p+- 


Pi 


<p{x',k'l) 


£+=o 


( 88 ) 


We isolated the term in brackets in the above equation because we need to specify how to handle the exponential at 
infinity mm- in the sense of principal-value distribution, we can write 


e ~ioop+ e~ iLp ' 

-x— = lim - 

p 7 L-hx 


x = 2?r ixd(pt) 

Pi 


(89) 


where x is a constant that takes different values, according to the prescription that we adopt for the denominator: 

(90) 


— 1 Retarded , 

X = { 0 Advanced , 

_ — \ Principal value . 


Again in the sense of principal-value distribution, instead, we have 

e -ioop+ = Q 


(91) 


which allows us to drop the first term in the bracket of ( 881 . This is not surprising, because that term comes from 
the —g pLI in the polarization sum, which appears in Feynman gauge also; then, this is consistent with the fact that 
the transverse gauge link does not contribute in the Feynman gauge. 

Now we can easily perform the integration over £ - , and then over x', k'_i_, to arrive at 


Ia = ~ 


*X e 


xP + 


(2tt) 2 [(1 - x)M 2 + k 2 ] 


S(p+) . 


(92) 
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As for term /g, which is now 


I B = -eP+ 


2 [ dt;-d 2 z ± 


l -loop' 

Jz _ 1J2' 


(27r) 3 J (2tt) 3 y/x'(1 — x') 


Pi 


■ x'ip{x ', k!1) 


«+=o 


(93) 


we notice that in any prescription (except for the advanced, where it is vanishing), it would become proportional to 
5(1 — x), thanks to (89); furthermore, the integration over would now bring a 5( 2 )(kjJ. Then, we can conclude 
that this term would be relevant only in the end point (x = l,kj_ = 0j_), and for this reason we will not consider it, 
consistently with what we did previously in the Feynman gauge. 

A similar argument can be applied to the term II. It can in turn be split into the sum of two contributions, of 
which the first one becomes relevant only at the end point, while the second one reads 


IIA = ~ 


eP+ 


dx' 


y’O^ki)<5(1 - x) . 


(2tt)3 J ~~ vV(l - *') pj 
In the sense of principal-value distribution, also in the above equation we can write 


(94) 


—— = 2 t n£5(p+) , 

Pi 

where in this case, according to the different prescriptions, the constant £ takes the values 

0 Retarded , 

£ = { 1 Advanced , 

I Principal value . 


(95) 


(96) 


We now need to fix a prescription in order to proceed further. In the retarded prescription IIa is vanishing, while I a 

xP + 5(1 — x) _ ie 2 


becomes, from (92) 


IA = 


1 


5(1 — x) . 


(2tt) 2 [(1 - x)M 2 + k 2 ] P+ ~ (2t r) 2 k 2 ± 

If we were to fix the advanced prescription, instead, we would find that Ia is zero, while IIa would read 


IIA = 


—ieP + 


dx' - 


5(1 — x') 


5(1 — x) = — - 


1 


[(1 - x') 2 M 2 + k 2 ± ] P+ v ' (27r) 2 k^ 

In the principal-value prescription, instead, the sum of Ia and IIa gives zero. 

We may summarize our resull 
in the different prescriptions, is 

Retarded: 


5(1 — x) . 


(97) 


(98) 


We may summarize our results by stating that the contribution of the transverse gauge link to the TMD /f (x, k^), 


1 


Advanced: — 


(27T) 2 k 2 
ie 


5(1 - x) , 

_L 

2 1 


(27t) 2 k 2 ± 
Principal value: 0 . 


5(1 - x) , 


(99) 


Adding the contributions in (99) to the results (78), (79) and (80) in the corresponding prescriptions, we can see 
that the transverse gauge link allows one to compensate for the extra terms that come from specific choices of a 
prescription for the regularization of the pole in the polarization sum, thus recovering the results obtained in the 
Feynman gauge. Similar results can be derived also for the other TMDs. We conclude that the transverse gauge link 
makes the evaluation of TMDs prescription independent in the light-cone gauge. 
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FIG. 7: Handbag diagrams related to the TMDs of the photon. 


VI. PHOTON TMDS 


We now focus on the derivation of the TMDs for the photon inside the dressed electron. We still consider the 
scattering of the electron off a lepton, with the situation this time described by the handbag diagram in Fig. [TJ We 
will now indicate with k = ( xP + , k~, k x ) the four-momentum carried by the photon. In analogy with gluon TMDs 
in the QCD case @210309, the relevant correlator between initial and final electron states for photon TMDs is 6 


J e ik <(P, S\F+*(0)F +i ^)\P, S) 


( 100 ) 


£+=0 


The held tensor in QED is already invariant under the gauge transformation ([5]); therefore in this case a 
Wilson line is not required. The situation is different in QCD, where the field tensor is not gauge invariant and gluon 
self-interactions are present. 

If we introduce the symmetry operator S, whose action on an operator O" is defined by 


SO" • - O" + O ji - Sij V O' 


( 101 ) 


we have the following decomposition of the correlator (1001 in terms of the leading-twist TMDs of a photon inside 


the dressed electron, which are obviously functions of (x, k^_) [42] : 
2 S j 




i4^\x,k ± ;S) =S z gJ L + 




IT » 


k, s 


m 


-L 7 
~9l T 


Cj Li h.3 

kr,S) = - + s s 


S k\ , S k\ 


VLFlVL ^J-7 , 
2m 2 1L + 


i Jk Q k 
± e ± d l 


2 m 


hi 


IT 


2 m 2 1T 


Sfcjei fc (2fefek x ■ S L - S K ± k 2 ± ) 
4m 3 




IT ’ 


( 102 ) 

(103) 

(104) 


where the antisymmetric tensor e" was defined in (14). For the photon TMDs we adopt the nomenclature of Ref. 


the letters /, g and h refer to the photon being unpolarized, circularly polarized or linearly polarized, respectively, 
while L and T still refer to the dressed electron being longitudinally or transversely polarized. Out of the eight 
photon TMDs implicitly defined in equations (102)-(104), four are T-even (fi,gJ L ,gJ T and h^ 1 ) and four T-odd 
C/rr , h i i. h iT and h i r )■ 


If we work in the liglit-cone gauge, the components of the field tensor entering the correlator (1001 reduce to 
F +l = d + A 1 . If we also switch to the light-front helicity basis and apply the LFWF overlap representation method, 


Indices i and j here refer to transverse Lorentz components. 
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we come up with 


$ 7M (j,kL,A) = 2( , 2 1 ^3 ^a,a; ((! -a^-kiJ^A,^ ((1 - z),-kj_) (fc)eA T ( fe ) ■ ( 105 ) 


Further development of ( 105[ ) yields 7 

t*](x,k ± ,A) = ly ( (1 2 -y> |2 6 


(A' 8a(P) I-y+ ^f7 + I (>1, + ™) ^-7* + jiff T + J “"( p ) ' 


Li 


(106) 


where the momentum of the inner electron is now p e = ((1 — x)P + The Feynman gauge result is completely 

equivalent, the difference in the polarization vectors being compensated by the extra diA + term in the field tensor. 

If one takes the proper combinations of helicities for the dressed electron and transverse indices, it is easy to recover 
from equations (102)-(104) the following results for the leading-twist photon TMDs: 


fi ( x , kj_) = 


9il( x >*±) = 


e 2 ky [l + (1 — a;) 2 ] + m 2 x 4 
(2tt) 3 x [k^ + m 2 x 2 Y 
e 2 kj(2 — x) + m 2 x 3 


(2t r) 3 

[k^_ + m 2 ir 2 ] 2 

2e 2 

2 x(l-x) 

(2tt) 3 

lib n 

[k^_ + m 2 x 2 ] 

2e 2 

-- —77 

2 (! - x) 


frr( x i k l) = k l) = h^{x,k\) = hJ T {x,k 2 ± ) = 0 . 


(107) 

(108) 

(109) 

( 110 ) 
( 111 ) 


Notice that the T-odd functions are vanishing: this is not surprising, since there are no contributions that can be 
interpreted as final-state interactions. Even if we go up to order a 2 , we find that T-odd TMDs still vanish. This is 
probably true to all orders, since it is related to the absence of a gauge link in the correlator. The vanishing of the 
Sivers function up to order a 2 is also consistent with the sum rule derived by Burkardt jT8l 1501 foil , since we found 
that the Sivers effect for the electron vanishes both at order a and at a 2 . 


VII. RESULTS 


In this section we discuss the results of our calculation for the electron and photon TMDs at 0(a ), which is 
strictly valid at (x,kj_) ^ (1,0). In Fig. § we present the unpolarized TMD /f (upper-left panel), and the TMD 
gl L (upper-right panel) for longitudinally polarized electrons in a longitudinally polarized electron, as functions of k 2 ± 
and for different values of x. In the lower panels of Fig. [8j we show the corresponding results for the combinations 
f +e = (/i + 9il)/2 and f~ e — (/i — 5 il)/2 , describing the probability to find the internal electron with spin aligned 


and antialigned to the spin of the parent electron, respectively. As discussed in Sec. IV A /f and gf L are given by 


different combinations of the squares of S- and P-wave components. In particular, /f is obtained from the sum of the 
contribution from the two partial waves, while in gf L one has the difference between the S- and P-wave contributions. 
Correspondingly, with the combination f +e we isolate the contribution from P waves to both ff and while f~ e 
gives the S-wave contribution which enters with a positive sign in and a negative sign in g\ L . We notice that for 
k' 2 ± -y 0, the P-wave contribution is vanishing, while the S-wave contribution reaches its maximum, with larger values 
at increasing x. On the other side, the falloff in kj_ of the S-wave contribution is very steep, and the contribution of 
the P wave takes over at increasing transverse momentum. This also means that, for vanishing transverse momentum, 
the configuration with the spin of the internal electron being antialigned with the spin of the target electron and 


7 


Notice that the denominator from the photon propagator does not need to be regularized, since we are not integrating over the photon 
momentum. 
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the spin of the photon being aligned with the spin of the target electron is favored, and at larger k\ the spin of the 
internal electron is most likely to flip in the direction of the spin of the parent electron. The spin-flip of the active 
parton at low transverse momentum is also responsible for the sign change in gf L , which occurs faster in k 2 ^ at larger 
values of x. 

In Fig. [ 9 ] we show the results for the transversity TMD h \, and the combinations h +e = (/f + /if)/2 and h~ e = 
(/■f — /if)/2, describing the probability of finding the internal electron with transverse spin aligned and antialigned 
to the transverse spin of the parent electron, respectively. Being chiral odd, /if involves a helicity-flip of the internal 
electron from the initial to the final state, which is compensated by a helicity flip of the parent electron in same 
direction. Since total helicity is the same in initial and final states, /if is diagonal in OAM, and, from the LFWF 


overlap representation in Eq. (67), we note that it is given in terms of the partial waves with L z = ±1 only, without 


contributions from the S-wave components. The P-waves enter in /if with different combinations with respect to the 
contribution in /f and gf L , as given by f +e . In particular, from the LFWF overlap representation in Eqs. (62), (63), 


and (67), we find that f +e contains the sum of the square of the L z = +1 and L z = — 1 components, while /if is given 


by the interference between the two P waves, i.e. 


f +e (x, kD = 


1 

2(2tt) j 


[|^,+i(^,k ± )| 5 


|^+_i(a;,k L )| 5 


( 112 ) 
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FIG. 8: TMDs /J(ai,k^) (upper-left panel), g^ L (x, kj_) (upper-right panel), along with their sum (lower-left panel) and their difference 
(lower-right panel). All distributions are rescaled by a factor 27r 2 /a, and shown as function of k^, and for different values of the longitudinal 
momentum fraction x: x = 0.1 (solid curves); x = 0.2 (dashed curves), and x = 0.5 (dashed-dotted curves). 
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FIG. 9: TMDs hf(x, k^) (upper panel), h +e (x,h 2 ± ) (lower left panel) and h~ e (x, k^_) (lower right panel), rescaled by a factor 27r 2 /a, as 
function of k^, and for different values of the longitudinal momentum fraction x: x = 0.1 (solid curves); x = 0.2 (dashed curves), and 
x = 0.5 (dashed-dotted curves). 


hi ( x i kj_) = 


(2tt) 


:Re 


*+i(x, kjjtf + 1 (z, kj_) 


(113) 


in Eqs. (36), (37) and (41), one also finds 


where we used the property (57) of the two-body component of the electron LFWFs. From the analytical expressions 


h\(3 


’ k D = YTx^ f+e ^ ±) 


2 \ 


(114) 


As a result, h\ has the same dependence on k\ as / +e , but its overall size is always smaller, particularly at smaller 
values of x. In the combination h +e , the positive contribution from the S-wave components of /® is evident at vanishing 
transverse momentum, while at higher values of kj_ the contributions of the P waves from /f and h\ become dominant. 
In the case of h~ e the P-wave contributions from /f and h\ partially cancel, especially at higher values of x , and we 
find that the low-fc^ behavior is dominated by the S-wave contribution of /f, while the P waves still control the tail 
at higher transverse momentum. As a consequence, at x = 0.5 the transverse spin of the internal electron tends to be 
aligned with the transverse spin of the parent electron in the full kj_ range, while for smaller values of x there is no 
marked preference between the parallel or antiparallel configuration of the two transverse spins. 


In Fig. 10 we show the TMD gf T for longitudinally polarized electrons in a transversely polarized electrons target 
in comparison with describing the momentum distribution for transversely polarized electron in a longitudinally 
polarized electron target. These distributions are characteristic effects due to transverse momentum, since they are 
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FIG. 10: TMDs g? T (x, k( ) (left panel), (x, k 2 ) (right panel), rescaled by a factor 2ir 2 /a, as function of k( , and for different values 
of the longitudinal momentum fraction x: x = 0.1 (solid curves); x = 0.2 (dashed curves), and x = 0.5 (dashed-dotted curves). 


the only ones which have no analog in the spin densities related to the GPDs in the impact parameter space J52] , and 
vanish after integration over k±. Since the LFWFs ’Pi and 4Q +1 vanish, the LFWF overlap representation in 
Eqs. (64) and (65) can be rewritten as 


9 it ( x > kj_) — 


2(2-7r) 3 k 3 


Re 


m 


(2^) 3 ki 


Re 


k R VP| _ x (x, kj.)^ + _ 1 (a:, kj_) 
k R *± >+1 (x, k_L)^+ +1 (x, kj.) 


(115) 

(116) 


We notice that gf T requires a helicity flip of the electron target that is not compensated by a change of the helicity of 
the internal electron, and vice versa h££ involves a helicity flip of the internal electron, but is diagonal in the helicity 
of the electron target. As a result, the two distributions require a transfer of OAM between the initial and final states: 
in the case of gf T , this comes from the interference of the S wave and the P wave with L z = +1, while for h££ it is 


driven from the S wave and the P wave with L z = — 1. The difference between the two P waves in Eqs. (49) and (50) 
leads to the following relations for the TMDs: 


9iAx,k 2 ± ) = -xh±l{x, ki 


(117) 


which explains the suppression of g^ T: in absolute value, with respect to hf£, especially at lower values of x. 

From Eqs. and omitting the contributions of , h£ e and hf£. which are vanishing at leading order in 

a, we can form the following densities for electrons of definite longitudinal or transverse polarization: 


p(x, k x , k y , (A, s_l), (A, SjJ) = ^ 


ft + AA gt L + ASj_^j.— 9it + Asyfcy— hj;£ + S ^ h\ 

m in 


(118) 


where we denote the transverse spin of the internal electron with sj_. 


In Fig. 11 we show the densities in the transverse momentum plane and at fixed x = 0.5 for a longitudinally 
polarized electron target and longitudinally polarized internal electron, with helicity in the same (/ +e ) or opposite 
( f~ e ) direction, along with the corresponding results for f f. As discussed above, the behavior of f +e is determined 
from the P waves, which rapidly increase starting from the dip in the center, reach a maximum at |fcjJ « 200 keV 
and smoothly decrease at larger values. On the other hand, the S-wave contribution to f~ e is maximum at the centre 
and rapidly falls off towards the periphery of the k x -k y plane. The interplay between the P-and S-wave contributions 
determines the pattern of /f = f +e + f~ e from the center to higher values of k±. In summary, the density of electrons 
in a dressed electron in momentum space, averaged over all polarizations, at x = 0.5, looks like a ring-shaped image, 
with a radius of about 200 keV. 

In Fig. 12 we show the transverse polarization densities h +e and h~ e . Qualitatively, they look similar to their 
longitudinal counterparts, but they are different in the details. The P-wave contribution is more pronounced in h +e , 
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FIG. 11: Density plots in the transverse momentum plane k x — k y and at fixed value of x = 0.5 for /+ e (upper left panel), f ~ e (upper right 
panel) and ff (lower panel), rescaled by a factor 2ir 2 /a. The legend in the bottom-right corner of each panel indicates the corresponding 
spin configurations: the grey and white empty discs refer to the unpolarized dressed and unpolarized bare electron, respectively; the circle 
and cross inside the discs stand for polarization along the longitudinal axes, in opposite directions. 


but the S wave is also present, and the density at kj_ = 0 is not zero. The S wave dominates in h~ e , but a P-wave 
contribution is also present, even though it is suppressed. Note that, in principle, the transverse polarization density 
could in principle be not cylindrically symmetric. However, the fact that and h£ e vanish in QED at this order 
makes the density symmetric. 

In Fig. [13] the distorting effect induced by the polarization is shown. In the left panel, we present the density 
in the transverse momentum plane and at fixed x = 0.5 for longitudinally polarized electrons in an electron target 
transversely polarized along the y direction. This is obtained by adding the dipole deformation due to the term gl T 
in Eq. (118) to the monopole distribution of /f. The term in g^ T features a significant dipole deformation along the 
direction of the spin of the parent electron, and is responsible for the sizable downward shift along y. An analogous 
observation can be made for transversely polarized electrons in a longitudinally polarized electron target, as displayed 
in the right panel of Fig. [l3j In this case, the effect of the distortion is even more pronounced and in the opposite 
direction, since the strength of the deformation is due to the term — h^£, with h^£ = — 2g < [ T at x = 0.5. 

We avoid showing plots of all the photon TMDs, but we take the chance to make a few observations. The functions 
/7 and gJ T are similar to their electron counterparts with the replacement x o (1 — x). The helicity distribution has 
no straightforward relation to the electron helicity distribution. At variance with the electron case, the combination 
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FIG. 12: Density plots in the transverse momentum plane k x — k y and at fixed value of x = 0.5 for h +e (left panel) and h~ e (right panel), 
rescaled by a factor 27r 2 /a. The legend in the bottom-right corner of each panel indicates the corresponding spin configurations: the grey 
and white discs refer to the dressed and bare electron, respectively; the arrows indicate polarization along the y direction. 
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FIG. 13: Density plots in the transverse momentum plane k x — k y and at fixed value of x = 0.5 for the sum (/f + g \ T )/2 (left 

panel) and (/^ + h^£)/2 (right panel), rescaled by a factor 27 t 2 /q;. The legend in the bottom-right corner of each panel indicates the 

corresponding spin configurations: the grey and white discs refer to the dressed and bare electron, respectively; the circle inside the discs 
stands for polarization along the longitudinal axes; the arrows indicate polarization along the y direction. 


/ _7 contains contributions only from P waves, while / +7 contains contributions from both P and S waves. The 
other nonvanishing photon TMD is the so-called Boer-Mulders function describing linearly polarized photons 
in an unpolarized electron msaES]. We plot it in Fig. [14] (left panel), where we observe that it is very large, 
especially at low x. The corresponding function for gluons in a proton has been the topic of a few articles in the 
last years [MUbfij . mainly because it does not require polarized targets and it can be generated through perturbative 


QCD corrections m , similarly to what happens in our QED calculation. In Fig. [Ml (right panel), we plot the 
combination (/7 + (k 2 — k 2 ) / (2m 2 )h^ 1 ) /2 corresponding to the density of photons with linear polarization along the 
y direction in an unpolarized dressed electron. The distribution has two peaks, shifted along the polarization direction 
by approximately ±200 keV at x = 0.5. 
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FIG. 14: Left panel: photon TMD h \ 7 (x. kj ). rescaled by a factor 2tt 2 /a, as function of k 2 , and for different values of the longitudinal 
momentum fraction x: x = 0.1 (solid curves); x = 0.2 (dashed curves), and x = 0.5 (dashed-dotted curves). Right panel: density plot in 
the transverse momentum plane k x — k y and at fixed value of' x = 0.5 for the combination (/(' + (fc 2 — /.: 2 ) /(2r» 2 )h { ~)/2, rescaled by a 
factor 27r 2 /a. The legend in the bottom-right corner of indicates the corresponding spin configurations: the grey empty disc refers to the 
unpolarized dressed electron, the white disc with arrow indicate a photon with linear polarization along the y direction. 


VIII. CONCLUSIONS 


In this article we analyzed the dressed electron at order a as a system effectively composed of an electron and a 
photon. In analogy with QCD parton distribution functions, we studied all twist-2 polarized transverse-momentum- 
dependent distribution functions (TMDs) of the electron and photon in the dressed electron, which describe the 
three-dimensional structure of the dressed electron in momentum space. 

We performed the calculation of electron TMDs in Feynman gauge and light-cone gauge up to order a and for 
i / 1, ki / fli. We explicitly checked the equivalence of the two gauges and the critical role played by Wilson 
lines in the gauge-invariant definition of TMDs. This result was already known for the QCD case, but the explicit 
QED calculation is particularly useful for illustration purposes. A nontrivial question concerns the contributions 
of transverse Wilson lines, which in our QED calculation at order a appear only at the end point x = 1. In this 
case, we showed how the choice of different prescriptions for the photon propagator is correctly compensated by the 
contribution of transverse Wilson lines. 

In light-cone gauge, we performed the calculation of TMDs in two different approaches, namely using Feynman dia¬ 
grams and using the overlap representation in terms of light-front wave functions (LFWFs), confirming the equivalence 
of the two approaches. 

We calculated photon TMDs using their gauge-invariant definition which, in the QED case, does not require the 
addition of Wilson lines. The results are in fact the same for Feynman gauge and light-cone gauge. 


The results of the electron TMDs are summarized in Eqs. (361 to (42). T-odd TMDs (h^ e and f^f) are zero, a 
result that holds true also at order a 2 due to the lack of contributing diagrams. Also the T-even TMD vanishes 
in QED at order a. 

The results of the photon TMDs are summarized in Eqs. (107) to (111). All T-odd TMDs (/{^, , hJ T ) 


are zero. We checked that this is true also at order a 2 , and we argued that it is probably true at any order a" 

In the last part of the article, we showed several plots that visualize TMDs in various ways, including multidi¬ 
mensional density plots of the structure of the dressed electron in momentum space. For instance, we showed that 
the density of electrons (and photons, by symmetry) in a dressed electron in momentum space, averaged over all 
polarizations, at x = 0.5, looks like a ring-shaped image, with a radius of about 200 keV. 

We observed that the inclusion of polarization causes several interesting effects. Different combinations of longitu¬ 
dinal polarizations make it possible to isolate contributions with different projections of orbital angular momentum. 
Some polarized distributions are not cylindrically symmetric. For instance, the density of linearly polarized photons 
inside the dressed electron (related to the TMD h^ 1 ) has a distribution with two peaks, shifted along the polarization 
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direction by approximately ±200 keV at x = 0.5. The shifts increase towards lower x. Transverse polarization of the 
dressed electron causes distortions in the density distributions of longitudinally polarized electrons, with shifts of the 
order of 200 keV at x = 0.5, decreasing towards lower x. 

Our analysis can be further extended by including the contributions from the end points and the treatment of 
infrared and rapidity divergences. Another possible direction of study would be to identify which observables are 
sensitive to the TMDs of a dressed electron. 
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